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Abstract: We study maps preserving the Heisenberg commutation relation ab — ba = 1. 
We find a one-parameter deformation of the standard realization of the above algebra in terms 
of a coordinate and its dual derivative. It involves a non-local "coordinate" operator while the 
dual "derivative" is just the Jackson finite-difference operator. Substitution of this realization into 
any differential operator involving x and -4= , results in an isospectral deformation of a continuous 
differential operator into a finite-difference one. We extend our results to the deformed Heisenberg 
algebra ab — qba =1. As an example of potential applications, various deformations of the Hahn 
polynomials are briefly discussed. 
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1 Introduction 

The Heisenberg algebra Q 



[a, b] = 1 (1) 



made its first appearance in physics, long before the birth of Quantum Mechanics, through its 
realization involving a continuous coordinate x and a dual derivative d/dx = d, the latter being 
the basic differential operator of analysis. General differential operators, in one dimension, are 
then expressed in terms of powers of d, multiplied by functions of a; — a wide class of physical 
problems leads to an eigenvalue equation for such operators. 

The reason underlying the predominance of this particular realization in physics is the con- 
tinuous nature of most spaces under study. Recently, however, there has been a growing interest 
in discretized versions of spacetimc or other, internal spaces. This sometimes originates in the 
need for numerical computation, as in, e.g., lattice QCD, but has also been proposed as a model 
of small scale structure. On another front, alternative realizations of (Q) have emerged in string 
theory (see, e.g., ||, Jf|). In either case, x, d ceases to be the realization of choice and, in several 
cases, discrete (finite difference) operators acquire preferred status. Maintaining the validity of (Q) 
makes the transition from the continuous to the discrete non-trivial. The need then arises for new 
realizations of the Heisenberg algebra in terms of discrete operators. Given such realizations, the 
differential operators mentioned above can be deformed by replacing the continuous realization by 
a discrete one — the non-trivial feature of such deformations is that they are isospectralQ. The 
process may be regarded as a quantum canonical transformation. 

There has been already a considerable amount of research in this direction (see, e.g., [if], 
| : |, |T| and refer ences therein), with the discrete derivative ds, defined by 

d s >f(x)=S- 1 (f(x + 5)-f(x)), (2) 

receiving traditionally most of the attention. It can be argued that this is due, in part, to the fact 
that the form of the canonically conjugate "coordinate" variable xg is known (see |T(| and Ex. 3.3) 



It is clear, from its definition, that ds can be restricted to the (equally spaced) points of a lattice. A 
second natural choice would be an exponential lattice, the corresponding finite difference operator 
being the Jackson derivative (or Jackson symbol, see Q), defined by 

d q >f(x)^^—x- 1 (f(x)^f(qx)). (3) 

The problem with this choice, and the main motivation for this work, is that the form of the 
canonically conjugate "coordinate" operator x q seems to be unknown. We solve this problem in 
Sec. U below while in Sec. |^ we study representations, in a very general setting. Sec. |] shows how to 
generate new canonical commutation relation preserving maps from known ones and Sec. ^ briefly 
extends the above results to the case of g-canonical commutation relations. In Sec. |^ two isospectral 
deformations of the Hahn operator are presented as a concrete application — the corresponding 
polynomial eigenfunctions are also supplied. 

2 The Jackson Derivative and its Canonical Conjugate 

Given the Jackson derivative d q , satisfying 

d q x - qxd q — I , —l<q<l. (4) 

iWe ocassionally refer to ([!]) in the sequel as the canonical commutation relation (CCR). 
2 We do not consider here the effect of boundary conditions. 



Canonical Commutation Relation Preserving Maps 



3 



One finds 
where 



d q >x n = {n}x n -\ {n} = — n = 0,1,2,... (5) 



d q >x n = d q x n \0) s ,a t . (6) 

The notation used in the r.h.s. of the above equation is as follows. |0) denotes the "vacuum" , a ket 
annihilated by derivatives, d q \0) = 0. The subscript of |0) is an instruction: express all variables to 
its left in terms of x, d q (already in this form, in this particular example), then use the commutation 
relation (Q) to bring the d q s to the right of the x's. There they are annihilated by |0) leaving a 
function of x only — this function serves to define the l.h.s., i.e., the action of d q on x n . For a 
general function f(x), defined as a Taylor series in x, the above relation leads to the alternative 
definition (|3|), which makes it clear that d q acts on the exponential lattice {x, qx, q 2 x, . . .}. d q can 
be realized as a pseudodiffcrential operator 

d q ~ -^—x-^l-q^^x-^A}, A = xd, (7) 

where d is the partial derivative w.r.t. x, dx — 1 + xd, also annihilating the vacuum, d\Q) = 0. For 
the reasons mentioned in the introduction, one would like to realize also an operator x q , such that 
d q x q — l+x q d q . Using the commutation relation Ax~ x = x~ 1 (A— 1) we find x~ x A n = (yl+l)" -1 ^, 
so that (q = e h ) 

D q = ^x-^l-q*) 

= l — x -^Y l l-A n 

l-o z - ; n\ 

n—l 



= -- — x-'Y, -(i + Ay^d 

1 — a ^— ; n\ 

n— 1 

= -r^—(l + A)- 1 (q 1 + A -l)d 

d q = -^—B- 1 {l-q B )d. i B = l + A. (8) 

It will prove convenient in what follows to use the notation [a;] = j^y, with [0] = 1. Notice that 
lim 9 ^i[a;] = 1. We rewrite <M) 



d q = {Bj- 1 d. 



(9) 



d q is of the form d q — f(B)d, f(B) = \BJ 1 . We look for x q in the form x q = xg(B). Then 

d q x q = f(B)dxg(B) 

= f(B)g(B)+xf(B + l)g(B + l)d. (10) 

The r.h.s. above should be equal to 1 + x q d q . One concludes that g(B) = f{B)~ 1 : i.e., 



x q = x[B\ = [Ajx 



(11) 
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The action of x q on monomials is 



x q >x n = [n+ ijx 



n+l 



(12) 



d q above acts on power series as a discrete derivative. We examine the corresponding interpretation 
of the action of x q . To this end, we introduce the Jackson integral operator S (see, e.g., B) given 
by 

S = {A}- 1 x. (13) 

Notice that {A} is invertible on the image of x, i.e., on x n , n = 1,2, . . .. Comparison with (Q) 
shows that d q S = 1 while (Sd q ) >x n — x n , n = 1, 2, . . . and (Sd q ) > 1 = 0. 

/(*) 




Using the expansion 

oo 

one finds 

oo 

S»f(x) = (l-q)J2<l n xf(q n x)' 



(14) 



(15) 



n=0 



q 3 x q 2 x qx 



Figure 1: The action of S on /, evalu- 
ated at x, gives the area under the dotted 
lines. 



i.e., S > f{x) gives the area under the dotted lines in 
Fig. |l| and converges to dx' f(x') in the limit q — > 1. 
Using the second of (pTj), we find 



(16) 



X q = xdS . 



In other words, the action of x q on f(x) consists in first 
producing the function f(x) = (dS)>f(x) and then mul- 
tiplying this latter by the classical coordinate x. 

Aside: We derive an alternative expression for f(x). In classical calculus one has (Rolle theorem) 

/(*) = (/>- + . (17) 



where {-) x denotes (classical) averaging in the interval [0, x], 
derivative w.r.t. x. Let now {■)% denote quantum averaging, 



(J)l 



-S>f(x) 



■ Jq dx' f(x') and f'(x) is the (classical) 

(18) 



Then 



(dS) > /(a 



8(l-q)'£q n xf( q n x) 

OO 

(l-9)E9"(/(^) + ^"/'(9" 



f{x) = {f)l + (xf)l , (19) 
which is the g-deformed ("quantum") analogue of (0). O 
Since d q , x q obey the CCR, one can define a quantum action > q , in complete analogy to the classical 



one, 



dq>q fiXg) = dgf(x g ) \0)x q ,d q 



(20) 
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where, in the r.h.s., d q is commuted past x q until it reaches the vacuum, where it gets annihilated 
- the remaining function of x q is, by definition, d q > f(x q ) (notice that the subscript of |0) instructs 
to express everything in terms of x q , d q ). This extends to arbitrary operators G(x q , d q ) acting on 
functions of x q , just like in the classical case. It follows trivially thatQ 

G(x q , d q ) > q f(x q ) = </> 9 (G(x, d) > f{x)) , (21) 

where <p q : x i— > x q , d i— > d q , is the q-deformation map. We note in passing that the quantum 
averaging operator, M q = ^S, is the inverse of B. Notice also that the product x q d q is constant 
in q (i.e., invariant under the deformation), 

x q d q = xd = A. (22) 

This implies that, when dealing with the g-deformation of a general differential operator, terms of 
the form x^d™, with m < n, can be expressed entirely in terms of A and d q , the action of which 
is simpler. The resulting q-deformed operator is a differential-difference operator. This occurs, 
for example, in the study of the hypergeometric operator. Notice also that the invariance of A 
permits the exponentiation of the infinitesimal generator of (j) q in a trivial manner. Finally, it is 
worth pointing out that the map tf> q admits a non-trivial classical limit, which preserves Poisson 
brackets. Indeed, with p, q satisfying {p, q} = 1, where {■, •} is the Poisson bracket, one can easily 
verify that {f(A)p, g/(A) -1 } = 1, where A = qp, giving rise to a wide class of classical canonical 
transformations that might in itself be worth exploring. 



3 CCR-preserving Maps and Adapted Bases 

Consider the Heisenberg-Weyl universal enveloping algebra H, generated by a, b, with [a, b] — 1. In 
the sequel we work with a certain completion TL of Ti, that allows us to deal with e.g., exponentials 
in a, b. From the discussion of the previous section (see (^l[)), we abstract a map cf> q : Tt — > "H that 
preserves the CCR[] 

<(> q : a^a q = {Bj^a, b^b q = b\B], (23) 

where B now is B — 1 + ba. All equations in the previous section depend only on x, d satisfying 
the CCR and are therefore valid for (x, d) i— > (6, a). Although we use the particular map (f> q given 
above as an example, we emphasize that our results below are general. 

For any pair (a, b) of abstract generators that satisfy the CCR, we say that the set n = 
0, 1, 2, . . .} is an adapted basis for (a, b) if a, b act on it as lowering and raising operators respectively 

a\n) = n\n - 1) , 6|n) = |n + l). (24) 



Example 3.1 A classical adapted basis 

One particular representation of the Hciscnberg algebra is supplied by the subalgcbra generated by 
{1,6}, an adapted basis being given by \n) — b n . The action of the Heisenberg algebra generators 
on an arbitrary power series f(b) is a > f = [a, f] and b> f = bf. □ 

Suppose now we are given a CCR-preserving map Q , where a denotes any parameters <j) might 
depend on, and we wish to find an adapted basis for the deformed generators it produces. A 

3 We have <f> q (G(x,d)j = G(x q ,dq) — the ordering of the x's and 9's in G(x,d) is immaterial, as long as the 
same ordering is used in G(x q , dq). 

4 We use the same symbol as for the map in (fell) 
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general solution to this problem is possible if we further impose the restriction that <p a be counit 
preserving (i.e., <p a (a)\0) = 0). It is worth emphasizing that this requirement, although rather 
natural, excludes nevertheless a number of familiar CCR-preserving maps, like the rotation from 
x, d to a* , a in the simple harmonic oscillator. Keeping this observation in mind, we proceed to 
the following statement: given any CCR and counit-preserving map <f> a , one can find in general 
an induced map <p a : \n) i— > \n) a that maps any adapted basis {\n}} for (a, b) to an adapted basis 
{|n) Q } for (a a , b a ). Indeed, to any function f(b a ) one can associate its b-projection f(b) given by^ 



f(b) = f(b a )\0) b , a . 



We now show that 




(25) 



(26) 



where / = f(b), G = G(b, a) and G a = (j> a {G{b, a)) = G(b a , a a ). We have 



G a >f(b) 



0) 



Ga/(6)|0>6,a 
G a f(b a )\0) b , a 
G a t> a f(b a )\0)b,a 

<\> a (G>/(6))|0) M 
(G>f)(b a )\0) bta 

GTf(b). 



(27) 



We comment briefly on the steps that lead to (|27|). The first equality follows from (g), taking into 
account that <f> a is counit preserving, so that we may put a in place of a a in the subscript of the 
vacuum. The second equality follows from (|25|). In the expression G a f(b a )\0) b , a we are instructed 
to express a a , b a in terms of a, b, and then bring the a's to the right etc... One can do this though 
in several ways. The one shown above involves first bringing all a Q 's to the right of the 6 Q 's, then 
substituting a a = a a (b,a) and bringing the a's to the right (this is equivalent to annihilating the 
a Q 's themselves, since <j) a preserves the counit). At this point one is left with a function of b a 
which is clearly G a t> a f{b a ). Finally, one substitutes b a = b a {b, a) and brings the a's to the right. 
Given an adapted basis {\n)} for (a, 6), we construct the set {|n) Q }, where^j 




and claim that it is an adapted basis for (a a , b a ). Indeed, 



Also, 



(28) — — 

o-a\n) a = a a > 6" = a > 6™ = nb n ~ x = n\n—l) a 



b a \n) a = b a bl\0) = bl +1 \0) = \n+l) a . 



(28) 

(29) 
(30) 



5 Notice that the b-projection of / depends on the particular deformation <fi a used — for simplicity of notation 
we do not show this dependence explicitly. 

6 The vacuum used in computing b n = b™\0) is the ket \n = 0). 
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Example 3.2 A quantum adapted basis 

Continuing our earlier classical example, we now turn to the realization of the Heisenberg algebra 
provided by the map <fi q . We take \n) — b n and find for (j> q (\n)) = \n) q 



\ll) q 



where [n]! = [1] [2] . . . [n] and [0]! = 1. 



= b n 

= « 

= (blB]T\0) 

= Ml\n), (31) 



Example 3.3 A second discrete realization 

Consider the pair of operators 

a s = Ma) = S-\e 5a -l), bg = Mb) = ber Sa . (32) 

One easily verifies that [as, b$] = 1 — this is the 5- realization of the CCR mentioned in the 
introduction. We take again \n) = b n and compute 4>s(\n)) = \n)s 

\n)s = 6J|0> 

=' be~ Sa be- Sa ...be~ Sa \0) 

= b(b -6) .....(6-( n - 1)6) e- n6a \0) 

= b(b-S)-...-(b-(n-l)5) 

EE b^. (33) 

For the falling d-factorial polynomials b s , defined by the next-to-last line above (also known as 
5-quasi-monomials 0), it holds 

n 

b { s n) = ^ s{n, k)S n - k b k , (34) 
fc=i 

where s(n, k) are the Stirling numbers of the first kind. □ 



Example 3.4 The q-exponential as a b-projection 
Consider the spectral problem 

d q > f(x) = \f(x) . 

Relying on (p6|), we look instead at the. equation dt> g(x) = Xg(x) and compute f(x) above from 
/ = g. We get g ~ e and, using x n — lnjx n , we find 

f(x) = g(x) 

= e Xx «\0) x , d 

= 

n! \n}\ 
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oo 



n=0 

= e q (\x) , (35) 

i.e., the standard g-deformed exponential is just e g (x) = e x . More generally, if 

a a = f-\B)a, b a =bf a (B), 

the eigenfunctions of a a are X^=o 6 ™' where f a (n)\ = f a (l)f a (2) ■ ■ ■ f a (n), f a (0)\ = 1. □ 

It is important to emphasize at this point the formal character of the above results. In particular, a 
given eigenfunction g(x) of some differential operator may converge for all x, while its 6-projection 
g{x) might have a finite (or even zero) radius of convergence as it happens, for example, with the 
(5-exponential 

°° 1 

e 5 (x) = e^ = Y d -ft ) (36) 

n=0 



that appears as the projection e Xx for the map <ps of Ex. 3.3 



4 Composition of CCR-preserving Maps 

If 4> a , <j)f3 are CCR-preserving maps then so is their composition <j) a o Considering only smooth 
maps, with a smooth inverse, one arrives at the notion of the group of CCR-preserving maps 
(CCR-PM). In the sequel we impose further the requirement that our maps preserve the counit. 
We can then use the fact that 

(j) a o fa — <j) a o fa (37) 
to compute the induced map of a composition of maps. We illustrate this in the following 
Example 4.1 Composition of <p q , <fi$ 

Consider the map <f> q discussed earlier. For each value of q, <fi q is an element of the CCR-PM. 
However, {4> q , — 1 < q < 1} is not a one-parameter subgroup since, in general, there is no q such 
that <j) qi o (j) q2 — <f) q . Notice also that [<p qi , 4> q2 \ ^ 0, for finite q\, qi. Similar remarks hold for <pg- 
For the composition cj> q o (f>g we find 

a q s 



(09 °M(a) 



= 5 



,sm~ 



b q s = {<t>q ° 4>s)(b) 
= Mbe' Sa ) 

= b(B}e- sm ~ la , (38) 



while for 4>$ o tfi q we get 



a.Sq = {4>8 O 4>q){o) 
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= <5- 1 [l + (5- 1 6(l-e-' 5a )]- 1 (e' 5a -l), 

b$ q = (<f> s o <j> q )(b) 
= MblB}) 

= be- 5a {l + S- l b(l - e~ Sa )} . (39) 
For the adapted bases that correspond to the above compositions we find 



\n)qS = <f>q O (f) S (\n)) 
n 

= ^ S (n,fc)^- fc 9 (6 fe ) 



fe=i 



Y,^n,k)S n - k lkjb k . (40) 



fe=i 



For example, 
Also, 



2 



\2) qS = 9 (b(b - 5)) = {2jb 2 -6b= —b 2 - 5b . (41) 



\n)s q = 4>8°4> q (\n)) 
= 0«([n]!|n)) 



[nf.br, (42) 



so that, for example, 

which should be compared with ([41|). □ 



\2) Sq = ^-b{b-5), (43) 



5 Quantum Canonical Conjugates 

Given x, y satisfying the g-Heisenberg algebra 

xy - qyx = 1 , (44) 

with —1 < q < 1. We say that x is the quantum canonical conjugate (QCC) of y (and vice versa). 
To complete our treatment of the map (j)s of Ex. |3.3| , we undertake here the determination of the 
QCC of ag. We work again with abstract operators a, b and remark that a q (given in (^3|)) and b 



satisfy (44), a q b — qba q = 1 . Notice that, the (classical) Heisenberg algebra admits the ^-involution 

b* = a , a* = b , (45) 

which we extend as complex conjugation to the parameter q, q* = q. Then, taking the * of ( ^4| ) 
(with a q expressed in terms of a, b, as in (||)), we find 

b*a q * - qa q *b* = 1 
- 1 ) - q (biB}- 1 ) a = 1. (46) 
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Up to now we disposed of the deforming map 4> q : a >— > a q , b > b, which, applied to a pair 
of operators satisfying the classical Heisenberg algebra, produces a pair satisfying the quantum 
Heisenberg algebra. Notice that it does so by leaving b invariant and only deforming a. What we 
have achieved in d46| ) , is to produce a second similar map (f>' q , which instead leaves a invariant and 
deforms only b: <\>' q : a l— ► a , b i— > a q * . We only need apply <fis to ( [ilf ) to get 

a s <f> s (a q *) - q<t>s(a g *) a s = 1 , (47) 

which identifies 4>s(a q *) = 4>s (6[_B] _1 ) as the QCC of as- 



6 Deformed Hahn Polynomials 

As an example of potential applications of our results, we present here various deformations of the 
Hahn operator and its eigenfunctions, the Hahn polynomials. We start with some definitions. The 
action of the Hahn operator Hs on a function fix) is given by 

H s >f(x) = <T 3 ( Ci 5 2 + c 2 5x + cix 2 ) f(x + 5) 

- <r 3 (c 4 <5 2 - 5{ Cl - 2c 2 + c 3 S)x + 2 Cl x 2 ) f(x) 

- S- 3 (S( Cl - c 2 + c 3 S)x - Cl x 2 ) f(x - S) , (48) 

where Cj, 8 are parameters. What distinguishes Hs is that it is the most general three-point finite- 
difference operator with infinitely many polynomial eigenfunctions p0| . The latter are called Hahn 

polynomials (of continuous argument) and we denote them by h^"'^^(x, N), where 

c 2 = N-2-/3, c 3 = -a-(3-l,c 4 = (l3 + l)(N-l), (49) 

and we have set, without loss of generality, 6 — 1 and ci = — 1. The corresponding eigenvalues are 

X k = 5- 1 c 1 k 2 +c 3 k, k = 0,1,2,... . (50) 

For particular values of their parameters and/or arguments, h^'^' S \x, N) reduce to the Meixner, 
Charlier, Tschebyschev, Krawtchouk or (discrete argument) Hahn polynomials (for details and 
their role in finite difference equations, see j|, ||, ||] and references therein). We will use the form 

h^ m {x,N)=j2l lX f, (51) 

where x^ is as in (|34|) and the 7, are known coefficients. 
It has been shown in jl0[ [ll| that Hs belongs to TC, 

Hs = ci(bsas) 2 {as + S^ 1 ) + c 2 bsa 2 + c 3 bsa s + c 4 a 5 , (52) 

where as, bs are given in (p2|), with (a, b) 1— > (d,x). We are now at a setting where our earlier 
results may be applied directly. First, we deform isospectrally Hs to H, by effecting the substitution 
(as,bs) 1— > (<9, x) in ( p2[ ) — a little bit of algebra gives 

H = d x 2 d 3 + [(ci + c 2 ) + a&- l x] xd 2 + [a + [cio- 1 + c 3 )x] d . (53) 



The cor responding polynomial eigenfunctions are obtained directly from (51), using the results of 
Ex. U, 

fc 

h^\x, N) (54) 



i=0 
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(the eigenvalues are, of course, still given by (Q)). A second isospectral deformation, involving 
the substitution (d, x) i— > (d q ,x q ) in (|53|), leads to 

H q = Cl A 2 {d q + S- 1 ) + c 2 Ad q + c 3 A + Ci d q , (55) 

[A = xd) with polynomial eigenfunctions 

fe 

h { ^ q \x,N)=Y.^- xl - ( 56 ) 

i=0 

Notice the ease with which h^"'^ (x, N) are computed, despite the highly non-trivial complexity 
of the differential-difference operator H q . Finally, we point out that the spectrum of H in ( |53| ) 
may also be g-deformed by effecting the substitution (d,x) i— » (d qi x) in d55| ) — one gets a finite- 
difference operator 

H q = a x 2 d q 3 + (ci + c 2 + ci^a;) x<9 g 2 + [c 4 + (ci<T 1 + c 3 )x] 9, , (57) 
with infinitely many polynomial eigenfunctions — the corresponding eigenvalues are 

X k = c l 8- x {k}({k-l} + l)+c 3 {k}, k = 0,1,2,.... (58) 

Note Added 

After this work was sent for publication, Professor C. Zachos, to whom we express our gratitude, 
pointed out to us that x q , d q can be related to x, d via a similarity transformation. Indeed, 
starting with the ansatz J7(A)" 1 5C/(A) = d q , it follows that dU(A) = d q U(A - 1) and, using (0), 
U(A) = IA}~ 1 U(A - 1), from which the formal expression U(A) = T q (A + l)/T(A + 1) can be 
derived. Here T q denotes the g-deformed gamma function, T q (x + 1) = {x}T q (x) (see, e.g., Q). 
Then, x q is computed as x q = U(A)~ 1 xlI(A) = U(A)^ 1 U(A — l)x = [AJx and the constancy of A 
in q follows trivially. 
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